The Origins of Phase Transitions in Small Systems 
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The identification and classification of phases in small systems, e.g. nuclei, social and financial 
networks, clusters, and biological systems, where the traditional definitions of phase transitions 
are not applicable, is important to obtain a deeper understanding of the phenomena observed 
in such systems. Within a simple statistical model we investigate the validity and applicability 
of different classification schemes for phase transtions in small systems. We show that the 
whole complex temperature plane contains necessary information in order to give a distinct 
classification. 
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PACS numbers: 05.70.F, 36.40.E, 64.60.C 



The thermodynamics of small systems, e.g. Bose- 
Einstein condensates in magneto-optical traps |jl|, ||, D, 
the nuclear liquid-gas transition observed by niultifrag- 
mentation in heavy ion reactions |4[ |^, g], and the 
solid- liquid phase transition of sodium clusters |Q, 
have gained increasing interest over the last few years. 
Because these systems are far away from the thermo- 
dynamic limit the standard tools for the description of 
phase transitions are not applicable and new concepts are 
needed. Within the last few years several classification 
schemes for phase transitions in finite systems have been 
proposed ||]. In this Letter we compare these classifica- 
tion schemes by means of a simple statistical models for 
atomic clusters and show that graveling transitions oc- 
curring in these models can only completely understood 
by considering the whole complex temperature plane. 

Among others Gross et al. have suggested a micro- 
canonical treatment |l^, |ll|, where phase transitions 
of different order are distinguished by the curvature of 
the entropy 5 = fee \n^l{E). According to their scheme a 
back-bending in the microcanonical caloric curve T{E) = 
l/dE^n{Q,{E)), i.e. the appearance of negative heat ca- 
pacities, is a mandatory criterion for a first order tran- 
sition. Caloric curves without back-bending, where the 
associated specific heat shows a hump, are classified as 
higher order transitions. 

We have proposed a classification scheme based on the 
distribution of zeros of the canonical partition function 
in the complex temperature plane . The classical par- 
tition function 
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can be factored into a product of the kinetic part and 
a product depending on the zeros Bk = Pk + i^k, with 
B-k = Bl of this integral function in the complex tem- 
perature plane, 
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It has been shown that phase transitions can be classified 
by a set of three parameters(a, 7, n), describing the dis- 
tribution of zeros close to the real axis, where 7 = tan u 
is the crossing angle between the real axis and the line 
of zeros, and a is determined from the approximated 
density of zeros 4'{t) ~ on this line. For infinite 
systems it has been exactly shown that a = 0, 7 = 
and Ti — Q corresponds to a first order phase transition, 
while a > corresponds to a higher order phase tran- 
sition For finite systems ri is always greater than 
zero reflecting the size of the system. The classification 
scheme can be extended to values of a < also being 
interpreted as first order phase transitions. This scheme 
sensitively reproduces the space dimension and particle 
number dependence of the transition order in Bose Ein- 
stein condensates JT4I and the first order nature of the 
nuclear multifragmentation phase transition [l5| . The dif- 
ferences between both schemes can be revealed within a 
simple statistical model for atomic clusters. The multiple 
normal-modes model describes structural transitions 
within small noble gas clusters by considering several iso- 
mers and the vibrational eigenfrequencies of the isomers. 
For a two isomer system the partition function can be 
written as 
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= exp(-/3£;i) -f p2 exp(-/3S2)) , 

where the uij are the normal modes of isomer i and the 
Ui are the permutational degeneracies of the isomers and 
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lie on a straight line and are equally spaced yielding 
a = 7 = thus implying a first order phase transition 
in any case (Ai? — E2 — Ei ) (see Figure 1). It is impor- 
tant to note, that with increasing system size the energy 
difference between the isomers will also increase, thus 
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FIG. 1: Logarithm of the canonical potential energy difTcrcnce expectation value log {{AE)) in the complex temperature plane 
for (a) P2/P1 = 50000, (b) P2/P1 = 5000, and (c) P2/P1 =0.5. The location of the zeros of the partition function are signaled 
by the sharp needles. In all cases the distributions of zeros indicate first order phase transitions. 
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FIG. 2: Microcanonical caloric curves for the multiple normal- 
modes model with energy difference between the isomers 
AE = 20. For a) A;'=10,12, and 15 and constant P2/P1 = 
50000 the back-bending is manifest for N = 10, can be tuned 
out by increasing the particle number, and disappears for A'^ 
as low as N =15. In b) for constant JV = 10 the back-bending 
can be tuned out by decreasing the ratio p2 /pi . 



Ti approaches zero. The microcanonical caloric curve 
T{E) = l/c^s ln(0(£')) for this model can be calculated 
via the inverse Laplace transform il{E) = {Z{/3)}. 
Fig. 2 shows that the back-bending arrogated in the 
Gross scheme for a first order phase transition can be 
tuned in an out by variation of the model parameters. 

The kinetic part of the partition function /3~(^^~^) 
plays the crucial role. If this is taken into account 
the micro-canonical caloric curves change dramatically, 
whereas this part has no effect on the distribution of zeros 
(the particle number dependence of the canonical parti- 
tion function is not only reflected by the kinetic part 
itself but also implicitly by the ground state energies). 
The change in the topology of the configuration space 



or equivalently configurations space regions with signifi- 
cantly chaniging vibrational entropies seems to be a nec- 
essary condition for phase transitions in small systems. 
Similary results have been pointed out by Pettini et al. 
Within classical statistical mechanics the kinetic part of 
canonical partition function is separable and the parti- 
tion function splits up into a kinetic and a potential part 
which can be handled independently. Within the micro- 
canonical ensemble structural phase transition might be 
washed out or hidden by the kinetic energy contributions 
to the entropy. 
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FIG. 3: Canonical specific heat reduced by the kinetic contri- 
bution for the same values of P2/P1 as in Figure 1. For values 
with p2 /pi larger than 1 the expected signals of a first order 
phase transition are seen. The value p2/pi=0.5 corresponds 
formally to a first order phase transition at negative temper- 
ature. This graveling transition exhibits at positive temper- 
ature a very weak hump in the specific heat (the graph is 
amplified by a factor of 25). 

A very interesting feature of the multiple normal- 
modes model occurs in the case where the isomer with 
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the lower ground state energy has a larger vibrational en- 
tropy (see Fig. 1 (c)). In this case formally a first order 
phase transition at negative temperatures occurs. The 
structural transition between the isomers, which occurs 
when the temperature is increased is accompanied by a 
drop in the vibrational entropy. This is a graveling transi- 
tion with a significantly smaller signal in the specific heat 
than that of the normal transition (see Fig. 3). Fig. 3 and 
Fig. 4 show i) that the zeros in the complex temperature 
plane sensitively detect phase transitions and ii) it is very 
important to use /3 as the natural variable since only this 
yields continuous pictures of thermodynamic properties. 
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FIG. 4: Specific heat in the complex temperature plane for 

P2/pi=0.5. The figure displays how the interplay of the zero 
and the pole of the specific heat influences the behaviour of 
the specific heat curve at positive temperatures. 

In conclusion we have found that the classification of 
phase transitions in small systems based on the curva- 



ture of the microcanonical caloric curves seems to be not 
rigorous enough to make distinct statements about the 
order. We have shown that the investigation of the whole 
complex temperature plane adds a valuable amount of in- 
formation in order to classify phase transitions in small 
systems. 
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